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Abstract. We demonstrate how to build a simulation of two dimensional physical theories describing topologically 
ordered systems whose excitations are in one to one correspondence with irreducible representations of a Hopf algebra, 
D(G), the quantum double of a finite group G. Our simulation uses a digital sequence of operations on a spin lattice 
to prepare a ground "vacuum" state and to create, braid and fuse anyonic excitations. The simulation works with or 
without the presence of a background Hamiltonian though only in the latter case is the system topologically protected. 
We describe a physical realization of a simulation of the simplest non-Abelian model, D(S3), using trapped neutral 
atoms in a two dimensional optical lattice and provide a sequence of steps to perform universal quantum computation 
with anyons. The use of ancillary spin degrees of freedom figures prominently in our construction and provides a novel 
technique to prepare and probe these systems. 
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1. Introduction 

Broadly speaking, topology deals with properties of spaces invariant under continuous transformations. Topology 
appears in a physical setting, typically, when a quantity X can be shown to take values in a discrete set {x\ ,xj, ■■■} 
when the parameters of a configuration space f vary over a continuous range, in which case J splits into sectors 
labelled by the discrete values of X, 

f = |J *i • (D 

a, ex 

these values providing a topological classification of configurations. Two configurations in different sectors cannot 
be transformed into one another continuously, in other words, quantity X is topologically stable and is, for instance, 
conserved during (smooth) time evolutions both in the classical and quantum worlds. Examples range from soliton 
theory in hydrodynamics, where soliton charges determine the stability of solitary waves, to the classification of 
instantons, the key to chiral symmetry breakdown in four-dimensional Yang-Mills theory. 

Topological phases of matter and lattice systems [ I ] are highly correlated phases whose order cannot be 
described in the local group symmetry and order parameter paradigm. Such systems exhibit a topological order 
usually associated with a gapped ground level with degeneracy dependent on the topological properties (typically, 
a Betti number) of the underlying spatial manifold — and, in two dimensions, with particle-like excitations with 
anyonic statistics, and gapless edge modes when boundaries are present. Fractional quantum Hall states exhibit 
such an order and provide an experimental setting where Abelian anyonic statistics have been shown to exist, and 
non-Abelian statistics are widely expected to appear for certain filling fractions. Topologically ordered lattice 
systems have also been theoretically constructed, and their experimental simulation complements their analytical 
and numerical study. 

Topological orders have recently attracted considerable interest in the field of quantum information, due to 
proposals to use topological phases of matter and lattice systems as quantum memories, whereby information is 
stored in topologically stable quantities, and quantum computers. In the most widely explored scenario, gates 
are performed on the stored information by creation, braiding and fusion of anyonic excitations. While it is not 
clear whether usual topological codes in two dimensions remain useful at finite temperatures or in the presence of 
noise [2], models in two dimensions have been proposed featuring string tensions between anyons [3], therefore 
confining anyonic defects — on the other hand, a kind of topological memories in four dimensions is also expected 
to show resilience to thermal effects. The anyon braiding paradigm of topological quantum computation hence 
deserves to be further studied and realised in the laboratory. 
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In this respect, it is instructive to consider, following DiVincenzo's example, the requirements an 
experimental setting must fulfil in order to implement anyonic topological quantum computation based on anyon 
braiding (see also [4]): 

• Existence: 

EO Existence of anyons. 

Anyonic statistics is used to act on code quantum states by manipulating anyonic objects. Quantum 

information may also be encoded in static configurations of anyons. 
El Controlled initialisation of anyons. 

Anyons may be present from the beginning in the system, or may need to be created; manipulation of 

anyons by transport or fusion may be required to bring the system to the desired initial configuration 

(e.g., a product state.) 
E2 Implementation of gates by anyon braiding. 

Braiding non-Abelian anyons results in the application of quantum gates. The statistics should be rich 

enough to provide a universal set of braiding quantum gates, in which case we speak of universal anyons. 
E3 Measurement of topological charge. 

This is the read-out part of a computation, and can be achieved in a number of ways, e.g., by using 

interferometry. Such methods may involve the controlled fusion of anyons. 

• Experimental feasibility: These are not independent of the existence requirements, since, e.g., the ability 
to perform gates imply that the life of an anyon is long enough to carry out the task, but their importance as 
criteria for practical implementations warrants a separate listing. 

50 Scalability of the implementation. 

A practical implementation should maximise the encoded information while ensuring control of the 
anyonic population. This leads to spatial efficiency as a design goal. 

51 Robustness of the implementation. 

For instance, the anyon lifetime should be larger than the gate application time, or than computation 
time for information-encoding anyons. 

The basic tools to perform anyonic TQC can be seen to be the controlled creation, transport, and fusion of anyons. 
The status of anyonic computation in different experimental settings is as follows: 

• Optical lattices: satisfy all of E0-E3 as argued in this article. 

• Quantum Hall effect: The existence of Abelian anyons as excitations of fractional quantum Hall systems is 
an established fact. Non-Abelian excitations have been argued to appear for filling fraction 5/2 and 12/5 — 
in the latter case, braiding would be universal. See recent experimental progress in [5]. 

• Josephson junction arrays: The theory was put forward in [6], leading to the experiments reported in [7]. 

• Photonic systems: See [8], [9] for recent experiments using entangled states of photons. 

Simulation of lattice systems using cold atoms and molecules in optical lattices offers an attractive setting 
where a high degree of control over the parameters of the system can be achieved. Minimal building blocks of 
topologically ordered systems have been considered in [10]. The implementation of Kitaev's honeycomb lattice 
model was discussed in [11] and [12]; braiding and interferometry experiments mediated by optical cavity modes 
were considered in [13]. Braiding in the Abelian phase was also considered in [14], however the degradation of 
anyonic visibility with this approach due to the perturbative treatment of the honeycomb Hamiltonian was pointed 
out in [ !], and a method to construct Kitaev's toric code from the cluster state was outlined in [16]. 

In this work we develop a practical, universal method, introduced in [ 7], to perform anyonic interferometry 
and arbitrary quantum computation tasks based on anyon braiding in spin lattices, using controlled operators 
applied by manipulating an extra species of ancillary particles to construct all needed primitives: anyon creation, 
transport and fusion. We illustrate the method using the D (53 ) quantum double model, whose excitations allow for 
universal quantum computation by braiding [IS]. Our construction is based on the spin lattice model introduced by 
Kitaev [ ] with the new feature that we give explicit protocols for manipulating anyonic excitations and read out 
of fusion products. We also describe a method to construct the relevant topologically ordered states efficiently in 
absence of a simulated topological Hamiltonian; however, while this construction is of interest in itself, we want to 
emphasise that the anyon manipulation procedure is independent and can be used whenever the topological phase 
is achieved, e.g., by simulation of the Hamiltonian. 
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2. Spin lattice models for the quantum double of a finite group 

Denote D(G) the quantum double of a finite group G = {gj}, which is a quasi-triangular Hopf algebra. An 
algebraic construction of these models was first given by Bais et al. [20] and for a brief introduction see 
Appendix A. We simulate a spin lattice Hamiltonian Hjo which is a sum of quasi-local operators and has localized 
particle like excitations are in one-to-one correspondence with the irreducible representations (irreps) of D(G). 
Consider a two-complex F, which is a cellulation of a two dimensional surface with vertex set V = {v,}, edge 
set £ = {e 7 }, and face set J = {//}. Particles with d = \G\ internal levels (qudits) are placed on the edges and 

physical states live in a Hilbert space 9i = H (c/)®l E where 9{{d) = C|0) H \-C\d—l). Particles on edges that 

meet at a vertex v all interact via a vertex operator A,,. Similarly, all particles on edges that are on the boundary of a 
face / interact via B f. We pick an orientation for each edge with e = [v ; , v&] denoting an edge with arrow pointing 
from vertex vj to Vk- The choice of edge orientations is not important as long as a consistent convention is used. 
We assume an orientable complex Y and each face / has an orientation (say, counterclockwise). The Hamiltonian 
is a sum of constraints chosen such that the ground states of Hjo are invariant under local gauge transformations 

T g (v)= J] L g (ej) J] R g Aej), (2) 

e;€[v,*] e;€[*,v] 

Here L g (e j),R g (e j) £ U(d) the d dimensional unitary group, are the permutation representations of the left and 
right action of multiplication by the group element g £ Gon the system particle located at edge e j. For the particle 
states we make the identification \ j) = \gj), where by convention |0) = \go) = \e), with e the identity element. The 
action of left and right group multiplication on the basis states is then L^j) = \hgj), Rf,\j) = \gjh). 
A suitable spin lattice model was provided by Kitaev [19]: 

ff T = -£A(v)-£i?(/) (3) 

/ 

where 

a(v) = ^I geG r g (v), 

Kf) = Eo^m^/I^^^I 

In the definition of B(f), the sum is taken over all products of group elements acting on a counterclockwise 
cycle of edges on the boundary of / such that the accumulated left action is the identity element e £ G (i.e. 
h(h(-\ . . .hih\ = e for the counterclockwise cycle starting at edge e\ and ending at edge eg). The function 
°f( e j) = ^1 according to whether the orientation of the edge is the same as(opposite to) the face orientation. 
By construction [A(v),A(v')] = [B(f),B(f')] = [A(v),B(f)} = 0. Furthermore, it is straightforward to verify that 
since A(v) is a symmetrized gauge transformation it is a projection as is B(f). The ground states of Hjo are then 

manifestly gauge invariant states. Excited states are described by violations of the local constraints A (v),B(f) and 

feel 

are particle-like corresponding to the irreps of D(G) labeled by Tv R ^ N ^ where [a] denotes a conjugacy class of 

G which labels the magnetic charge, and R(N^) denotes a unitary irrep R of the centralizer of an element in the 
conjugacy class [a] which labels the electric charge. Note there is an arbitrariness in how one picks the fiducial 
element of the conjugacy class. However, N ghg -\ = gNhg~ l so that the centralizers for the elements in a given 
conjugacy class are isomorphic (they are equal up to a gauge transformation), and we can index them just by the 
conjugacy class. 

We focus on a specific two-complex F which is a square lattice with boundary. For any two complex with 
boundary and without holes, there exists a ground state \GS) such that Hto|G5') = — | + \f DUGS') and it is 
unique (for the argument see e.g. [21]). The convention for edge and face orientations is shown in Fig. 1. We 
slightly abuse notation by labeling the particles according to location relative to a face index fj and vertex index 
Vj (see Fig. la). For instance, a vertex ancillary particle at vertex vy will be labeled vy and a face ancillary 
particle at face /y will be labeled /y. The system particle on edge e = [vy , /] will be labeled e/. j; i./. When we 
are referring to the actual spatial locations / and v it will be made clear. 

In Table 1 we give an algorithmic procedure to prepare the ground state of the Hamiltonian Hxo for an 
arbitrary finite group G. We begin with all system particles and face ancillae in state \e) = |0). This guarantees 
that the initial system state satisfies the zero flux condition, i.e. B(f)\\\t)s = |\|/)sV/. All vertex ancillae are 

prepared in the state |0) where \J) = -4=l[=o 1 e 2llijk/lGl \k). In the case G = la, this algorithm produces the 

Vl G 

ground state of the planar version of Kitaev's toric code [19]. For that model all operations can be done with 
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By single particle measurement, prepare the initial state I^V) = ®„ eE |0) e <8> v6 y [0} v (8) f e jr |0) /. 
for k = : m — 1 
for j = : n 

Apply the unitary W(vj^) where W(v) = £/, e c Wvv^l ® ?a(v). 

Measure ancilla v/j in the basis {| j)} (probabilty for outcome j is 1/|G|). 

For the outcome \j) apply the single qudit operation Zi(ej fcjfi+l ) where Z-'(e) = L[=o ' e' 2,t ^/l G l |fc) e (fc|. 

j++ 

end 

k++ 
end 
j=0 

for j = : n — 1 
Apply the unitary VK(vj.i). 

Measure ancilla v,-£ in the basis {\j}} (probabilty for outcome j is 1/|G|). 
For the outcome \j) apply Z j (e jtm; j + ^ m ). 

j++ 

end 

Table 1, Algorithm: Ground state-synthesis. An algorithm for preparation of the ground state | GS) of Hjo over a finite 
group G on an (n + 1) x (m+ 1) square lattice with boundary satisfying B(f)\GS) = A(v)\GS) = |GS)Vv £ t^,/ £ y. 
The algorithm works by beginning in a +1 co-eigenstate of all face operators B(f) and applying vertex ancilla 
assisted projections A(v). Using single particle operations conditioned on measurement outcomes of vertex ancilla, 

this algorithm outputs the state ll'out) = rLei'Vv [ / . Y,hea Ti,(v)] fin)- No projection A v is needed because for 

\ \J \G\ 

a two complex T with boundary, IIvev^Mk) = nvey\v„,„^( v )k) 8 ' E '- Each controlled gate operation requires 
0(\G\) elementary single and two qudit operations and the algorithm has complexity 0(\G\nm). Since [Lj,,R^/] = 0, 
the operators performing local gauge transformations commute: [W(i>), W(v')] = 0. Up to the last column then, all 
columnwise unitary operations W(v) and subsequent correction gates Z'(e) can be done in parallel. For the last 
column the correction gates do not commute with the operators W(v) hence the operations are done serially so that 
the computational depth of the algorithm is 0(\G\(m + n)). 

qubits, and the permutation rep of the group is L e = R e = I2 and L gl = R gl = (f. Controlled operations involve 
only CNOT gates and the correction gates Z' = (a*-)- 1 . Another scheme for constructing this ground state using 
single qubit measurements and feedforward on a prepared cluster state is given in [16]. 

Our algorithm has a computational depth of 0(\G\(m + n)) but one might wonder if a faster ground state 
preparation procedure is possible. The answer is no if the initial state is uncorrected and the available set of 
operations is quasi-local. The reason is that the final state has global correlations that are created by quasi-local 
operations. In our algorithm these operations are measurements but they could also be adiabatic turn on of the 
summands of Hjo- The time scale to perform the quasi-local operations (here the measurement time of A(v) 
establishes a light cone for the flow of correlations. In [22] it was shown by an application of the Lieb-Robinson 
bound that the minimal time to prepare a topologically ordered state beginning in a completely uncorrelated state 
is of the order of the length of the correlations. Since the correlation length scales as the linear dimension of our 
lattice, our algorithm is essentially optimal. 

3. Simulation ofD(5 3 ) 

3.1. Ground state preparation 

The algorithm Ground state-synthesis can be built using the tools for single- and two-qudit gates described above. 
The controlled gauge transformations W(v) are built up by applying L/, or R h -\ to left,right, top, or bottom 
neighbors of the vertex v controlled on v begin in state \h). Notice that unitaries L/, and R/, are maximally sparse 
in the logical basis, hence using the QR decomposition for simulating unitaries [23], each such operation can be 
built using a small number two qudit diagonal phase phase gates U = g'^Mw'MW-Bwi conjugated by single 
qudit Givens rotations. $. Measurement of vertex ancilla are performed in the Fourier basis \J) V , hence a one 
qudit Fourier transform F v is needed before measuring in the logical basis. That operator is not sparse but can be 
constructed with fewer than 13 steps using parallel Givens rotations [23]. 

% Using the software package in [24] we find a construction for the operator Y.hes^ \h)(h\ ® or Y.hes^ \h)(h\® Rh in 37 controlled phase 
gates. This count may be reduced to 8 controlled phase gates using a qutrit-qubit encoding for the group elements (see Appendix B) 
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Figure 1. A spin lattice model for topological order, (a) The system particles (blue) reside on the edges of a two 
complex r. In a Hamiltonian formulation, all particles on edges that meet at a vertex v interact as do all edges that 
surround a face /. The orientation of edges and faces is indicated. An ancillary set of particles (red) is placed on 
the vertices of the of the complex T and the faces of T (equivalently, the vertices of the dual complex f ). The ancillary 
particles on the vertices(faces) afford a handle on operations with system particles to create and guide electric(magnetic) 
charges which are indicated by diamonds(squares). The braiding of one member of a vacuum magnetic charge [£] pair 
around one member of a vacuum electric charge R pair is shown: ?R_ |0m;(v3 1 i,/3 l i),(vi i i,/i i i)}|l*;(v4 l 4,Vi i 5)) — > 

7 =E< £ Ml( i,sf . < );^,i./3,i)>l( pSf .«" , );(»i,i./i,i)>WfliK4,vi,5)). 



3.2. The particle spectrum 

The 8 irreps for D(Ss) are listed below with their corresponding quantum dimensions: 
IT^+ d = 1 (vacuum) 
IlL^ , n|j d = 2,3 (pure magnetic charges) 
Tv^_ , nj^ d = 1 , 2 (pure electric charges) 
Ilpj , lip! ' d = 2,2,3 (dyonic combinations) 

A complete derivation of the fusion rules for this model is given in [26]. 

Magnetic charges are labeled by conjugacy classes of the the group G. Recall the conjugacy class C/, is 
defined C/, = {ghg~ \g € G}. For S3 we label [e] =C e = S3, [c] = C c+ , [t] = C tQ . The magnetic flux across a face is 
given by the ordered product of group elements represented by the basis states of the edges surrounding the face. 
The order is taken along a closed counterclockwise cycle, beginning at an origin v. Except for the trivial flux e 
case, the origin from which the product is taken is also important as not labeling the origin is equivalent to only 
specifying the conjugacy class or magnetic charge. The projector onto states with magnetic flux I at face / as 
computed taking a connected counterclockwise cycle around / with the origin at vertex v is: 

Be(v,f)= I ®e^f\h7 f{ % k {h k ° M \ (6) 

{IL lt e3/A*=^ko=[v,*]} 

The projector onto states with magnetic charge [I] at face / is then the sum over all fluxes in the same conjugacy 
class: Y*ee[e]B(i(v,f), in which case the specification of origin is unnecessary. The identity element e is its own 
conjugacy class, hence we use the convention B e (v,f) = B(f). 

Electric charges are labeled by irreps of the centralizer of a conjugacy class of G. The centralizer is defined 
N h = {g<E G\gh = hg}. For S3, we have N [e] = S 3 ,N [t ] = {e,t } = Z 2 ,N [c ] = {e,c+,c_} = Z 3 . 
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For a finite group G every representation is equivalent to a unitary representation. Consider projection 
operators onto subspaces belonging to the unitary irreducible representation, or unirep, R with matrix 
representation/? and dimension \R\: 



\R\ 



(7) 



i'GG 



This projection operator satisfies: 

^vP^ = ^R.R'^.^x- 
as is verified by taking the inner product over G, 



P R P R' 
r p.v r K,X 



\R\\R'\ 
~W 

\R\\R'\ 

~W 

\R\\R'\ 
~W 

\R\\R'\ 



La'gg Ly -§f 5 k,v 5y i;U [R(ti)*] rX ti 



= j§S^S K>v E V6G [i?(/0*Uft' 

= &R,R'&k,vP^x 

We can obtain a set of |/?| orthonormal basis states \Rx) which satisfy 
rf y \Rx)=8R,R,8 v .x\R fJ ) 

by applying the set P R V over all /j for v fixed onto some vector \%) with P R v \x) ^ 0, and normalizing. 

For pure electric charges in D(53), we are interested in unireps of N^ = S3. There are three unireps: the one 
dimensional identity rep R^, the one dimensional signed rep R^ , and the two dimensional rep R^. The projection 
operators are 



1 



(e+to+h + t 2 +c + +c_) 



for the identity irrep, and 



r- 1 , 
o 



+ c + +c_) 



(8) 



(9) 



for the signed irrep, and: 

pR 2 








1 

1 

% 

^ 



'0- 



(10) 



for the 2 dimensional irrep 7? 2 with ^ = e' 2,t / 3 . In the context of the lattice spin model, charge at a vertex v 
will correspond to applying the components of these projection operators, with group action being local gauge 
transformations at vertex v, to some system state \ X V) 



P R t(v)\V)/y 

Pg6(v)W/\ 
< 2 (v)|^>/\ 
P?t(v)\V)/ x 

P??(v)\V)/x 



j. [T e (v) - T tQ (v) - T tl (v) - T h (v) + T c _ (v) + T c+ (v) 



-^[r e (v)+^r c+ (v)+^ c _(v)]|¥) 
^[r <b (v)+^r fl (v)+^%(v)]|'J0 

^[T l0 (v)+^T tl (v)+^T t2 (vm 

j s [Te(y)HTc + (y)H"Tc(ym 



(ID 



We obtain two copies of a basis for the 2 dimensional irrep R2, namely {-foo^io} an d the conjugate basis 
{^oo = Pii'Pio = ^oi 2 }- The system state |*P) here can be interpreted as the state at some intermediate stage 



of the state synthesis algorithm before a gauge symmetrization at vertex v has been performed. 
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For the dyonic combination with flux [c] the charges are labeled by unireps of the N[ c ] = Z3. There are three 
unireps of Z3: R®,R\, and R\, and since it is an Abelian group, they are all one dimensional. The projection 
operators for N Cp are 

P$ = |(e + c p +c_p), 

P R \ = \{e + %c p + %*c- p ), (12) 
P R \ = i( e + ^c p + ^_p) 

For the dyonic combination with flux [t] the charges are labeled by unireps of Afy = Z2. There are two one 
dimensional unireps of Z2: R\ and R\. The projection operators for N tj are 



P R 1 = \{e+ti), 
P R 1 = \{e- ti ), 



* f"T (13) 



Notice that there is a relation between projection operators for the centralizer of the identity and the products 

P^M'ocforoce [a]: 

P«!i(c + + c_) = cos©^, 
P R \\{c + + c-) = cos©J$ 

pR\ l/*._L»._L*.l — vR 



(to+h + tz) = P h 



P R *\(t» + h+t 2 ) = -P R T, 

In the spin lattice model electric(magnetic) charges correspond to violations of the local vertex(face) 
constraints and the state of a dyonic particle at vertex and face location (v,/) will be denoted: 

|(^ v ,g);(v,/)) 

where P R V labels an irrep of G as above and g is the flux at face /, evaluated taking a counterclockwise cycle with 
base point at v. For example, the state of magnetic flux I particle located at face / and its anti-particle located at 
face /' (with magnetic flux evaluated with respect to the origin v) is 

\{p R +,£y,(vj))\(p R +,r 1 y,(v,f)). 

For our model most excitations created in the bulk of the lattice appear as particle anti-particle pairs (such 
that the total charge of the pair is zero). Although single particle excitations can be made by creating them at the 
boundary. We will describe digital simulations of braiding of pure charges and pure fluxes. 



3.3. Anyonic dynamics 

Before deriving a sequence of operations to create and move anyons in the spin lattice let's review the rules for 
braiding charges and fluxes in anyonic models. We will write \a) to represent a magnetic flux of value a and 
|a,fl _1 ) for a flux anti-flux pair. For electric charge pairs we have one charge that transforms under the irrep 
R and the anti-charge which transforms under the complex conjugate representation R*. We introduce the bases 

{ «}jfJo 1 , { |v) r* }^! 1 on which the representations act and write a generic state of an electric charge anti-charge 
pair as a \R\ x \R\ matrix: 

\M R ) = -^=Y<Ky\v)R®\v}R* 
V\ K \ <".v 

with the normalization chosen such that ^ v |M^ V | 2 = \R\. 

Interchanging two fluxes, a left flux a and a right flux b in a counterclockwise sense is described by the action 
of the monodromy operator 3^ : 

H.\a)\b) = a\a}\aba^) = \aba~ l )\a) 

where a is the particle interchange operator. Squaring the monodromy operator gives the action of braiding two 
fluxes 

1l 2 \a)\b) = Kab^iaby^laba- 1 ) (14) 
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Braiding a flux b around a flux anti-flux pair (a, a ) is equivalent to braiding first around one then around the 
other (we can order the particles left to right (1,2,3)) 

^A^ul^hO = \b)\bab-\ba- x b- x ) (15) 

If \bab~ 1 ,ba~ l b~ l ) ^ then we say that \a,a~ l ) has magnetic charge. For each conjugacy class [£}, there 

is one unique chargeless state defined: 

l°M) = ill^- 1 > (16) 

v ll^JI te[(\ 

Electric charges moving past each other have no effect, only the braiding of fluxes around charges has an 
effect. Specifically, if we braid a flux \h) around one electric charge in the pair \M R ) we obtain: 

Hl 2 \h)\M R ) = \h)\R{h)M R ) (17) 
and if we braid around the anti-charge we obtain: 

H^\h)\M R ) = \h)\M R R{h- 1 )) (18) 
where R(h) is the matrix representation R of the group element h. Braiding around both charges is a conjugation 

^li®^h)\M R ) = \h)\R(h)M R R(h- 1 )) 
For each irrep R there is one unique fluxless state that is invariant under conjugation: 

|1|«|> = -J=2>)s®|aV. (19) 

3.3.1. Magnetic charges in T)(S^) Consider the creation of the chargeless magnetic flux pair for the conjugacy 
class [£] 

\Ow;(y, +1J+1 ,fij),(v i+1 j +l ,fij + i)) = -4= L |(^^r 1 );(v mj+1) /;- J ))|(A^);(v /+1J+1 ,/;- J+ i))(20) 

V li* J I te[i\ 

This state can be created by starting in the ground state \GS) and acting on one edge which is a shared boundary 
of the two faces: 

\0[ey,(vi+i,j+ufi,j),(vi+i,j+ufij+i)) = mm 52 R (( e iJ+i;i+iJ+i)\ GS ) 

v\m\ee[e] 

Note that the right multiplication operator 7?^ (e, ; commutes with all vertex operators except A (v/+i 

However, the sum of R( over all elements of the conjugacy class does commute with it. To reduce clutter here, we 
write v = V;+ij+i and e = e/.j. i./. i- then: 

^ m A{v)R c {e)A{v) = ^^fsWfj-W^e^VW^WV'W] 

= \^I, g , g >f g {v)t g ,{v)® \Z m R g -,{e)R g - Hg {e)R^{e)) 

= ] ^L g ^T g (v)T g iv)^\l ie[ei Re(e)R g - l (e)R gl - l (e)} 

= j^iL(e[C]Re(e)L g . g 'T s (v)T g/ (v) 

where T g (v)R g -i(e) = T g (v), since A(v)\GS) = \GS). Therefore, only the face constraints B(fij),B(fjj + i) are 
violated. Starting from the vacuum (ground) state \GS), this state is created by preparing the ancilla fij in the 
state |0 M )/y where |0 M ) = -4= Zee{(\ l^>. applying the two qudit unitary 

and measuring the face ancilla in the basis {\k [(] ) =Z* m \0 [(] )} l ^ Q \ where =^ me[ t ] e i2llkm l\^\\l m )(t m \ (where 
we have labeled the group elements in [I] = {/q, ... ,^]|-i})- For the outcome Om the target magnetic charge state 
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is created. Otherwise for outcome km, we need a correction step. To do this prepare the ancilla fij in the state 
\e)fij. Apply the controlled operation A(vi + \j + \,fij) where 

A(v,/)=£B,(v,/)®L,(/). 

g&G 

which maps the ancilla / to state \g) f when the flux at / (with respect to the base point v) is g. Such a 
controlled operation can be decomposed into elementary two qudit controlled rotation operators with each edge 
ek surrounding / as a control and the ancilla as the target, viz. 

A(vj + i .j+i ,fij) = A^ ght (fij)A bottom (fij)Ai efl (fij)A to p(fij) 

A(v i+ 1 j, fij ) = A top ifij) Arfght (fij ) Abottom (/i'j)Aleft (fij ) 
M v iJ , fij ) = A left (fij ) A to p (fij ) Arfght (fij ) Abottom (A./ ) 

A(V; J+ 1 , /j'j' ) = Abottom (fij ) Aleft (/ij ) A top (fij ) Arfght (fij ) 

where 

Aright (A;) = Eg£G l^esy+i^ij+i (g| ® V (/*»>)]> 

A bo ttom(/i,;) = \LgeG \g)e, r , J+l (g\ ® L g -l (fij)} 

Aleft(Ay) = lLgeG\g}e tJ ; MJ (g\®L g (f iJ )], 

Atop(/ij) = Eg6Gk)«W-ij i H.ij + i(*l®V 1 ^w)]- 

The convention here is that one takes the product of controlled left multiplication operations on the face ancilla 
along a counterclockwise path around / beginning at vertex v where the operator applied to / is L gT t depending 
on the orientation of the edge relative to the face. Next apply the single qudit phase gate [1 | G | _ | ^ | © Zp | ] (fij) and 

finally apply A(fij) to disentangle the ancilla from the system. 

We now describe in detail how to move magnetic charges from one face / to an adjacent face /'. Essentially 
it involves coherently mapping the value of flux at / to the face ancilla / and applying a controlled operation on 
an edge e G dfUdf' (the shared boundary of the faces /,/"')■ After this controlled operation the face ancilla / is 
disentangled from the system by mapping the flux at face /' to ancilla /' and performing a controlled operation 
between ancillae /,/' and finally reversing the mapping on /'. In this protocol we are careful to demand only 
single qudit and nearest neighbor two qudit interactions. 

Here is the protocol to move a magnetic flux one face unit to the right: 

| ,x)\ (v i+ 1 ,;+ 1 , ftj)) -> | (J*i , x) ; (v /+ 1 j+ 1 , fij+ 1 ) ) . 

• Prepare the face ancillae fij,fij+\ in state \e) 

• Coherently map the state of the flux at face fij to the ancilla fj via A(v,-+i j+\ ,fu) 

• Apply the controlled unitary Y(f ;,ei ;+i;j+i j+i) where we define 

Y (fA = I ZheG\h)f(h\®R h (e) o f (e) = +\ 
U, ) \ LheG\ h )f( h \® L h- l ( e ) o f (e) = -\ 
The remaining steps disentangle the ancilla from the system. 

• Map the flux value at the face fj+i to the face ancilla fj+i by applying A(vj+ij+i,/jj+i) 

• Swapquditse, J+1;i+1J+1 and/ u+1 

• Apply the unitary u = £ g€G \g) fl J+l (g\ ®L g -i (fij) 

• Swap qudits e, J+ i ;i+ i j+i and/y+i 

• Apply A(v i+ ij +u fij + i)- 1 . 

This entire process respects superpositions over flux states and can therefore be used to propagate magnetic charges 
around the lattice. A simple adaptation allows a magnetic flux move one face unit to the left: 

\(p R t, x y,(vij,fij)) - |(P 8 f,*);(vy J / 1 j_ 1 )). 

• Prepare the face ancillae fij,fi-ij-i in state \e) 

• Coherently map the state of the flux at face fij to the ancilla fij via K(yu,fij) 

• Apply the controlled unitary Y(fij,eij- t i+\j). 

• Map the flux value at the face fij-i to the face ancilla fu-\ by applying A(vij,fij-i) 



CONTENTS 



11 



• Swap qudits c U:i . u and fjj-i 

• Apply the unitary u = £ gGG \g) fi] _ l (g\ ®L g -i (fij) 

• Swap qudits eu-j+ij and fij-\ 

• Apply AfajJij-i)- 1 . 

Similarly, to move a magnetic flux one face unit down: 

|(J^,*);(v y+ i,/y)) - l(^,*);(vu+i»/i-u)>- 

• Prepare the face ancillae fij,fi-ij in state \e) 

• Coherently map the state of the flux at face fij to the ancilla fij via A(v; j+\ ,fi ,•) 

• Apply the controlled unitary Y(fi j , eij-u+i). 

• Map the flux value at the face fi-ij to the face ancilla / by applying A(v,- ,•) 

• Swap qudits e ( j ; y+i and 

• Apply the unitary u = £ geG (g)/^, ®L g -\ (fij) 

• Swap qudits e, j;i;/ ' + i and fi-ij 

• Apply A(vy+i,//-i,y) -1 . 

Finally, to move a magnetic flux one face unit up: 

x);(v ;+u , fij)) -» |(/^ jJC );(vh-u,/j + ij)). 

• Prepare the face ancillae fu,fi+i ; in state |e) 

• Coherently map the state of the flux at face /jj to the ancilla /; y via A(v,-+i j,fij) 

• Apply the controlled unitary Y (fij , e,+ 1 1 j+ 1 ) . 

• Map the flux value at the face fi+ij to the face ancilla fi+ij by applying A(v;+ij,/;+i j) 

• Swap qudits e (+ ij ;!+ ij+i mdfi+ij 

• Apply the unitary u = £ g€G \g)f i+1J {g\ ®L g -i (fij) 

• Swap qudits e (+ ij; !+ ij+i and fi +lJ 

• Apply A^i+ijJi+ij)' 1 . 

Fusion of a magnetic charge particle anti-particle pair can be measured by using controlled operations to bring 
the constituent charges in conjugacy class [£] adjacent to one another at faces (/,/') with shared edge e and 
applying first the operator A(v,/) = L«gg^( v >/) ® L g (f) followed by the operator Y*geG \s)f{s\ ®L g -i (e) then 
measurement of the ancilla / in the basis \gyr\) f. The probability to obtain the outcome O^j equals the probability 
for the pair to fuse into the vacuum. 

3.3.2. Electric charges in D(Si) Consider the following electric charge R(N([e]) particle at vy and its anti- 
particle pair at \'ij+i 

\(^e);( Vi j ,-))\(^,e);(v iJ+1 ,-)), (21) 

which is a labeling of the state by basis states of the irrep R and it's conjugate R*. The index p labels different 
copies of higher dimensional irreps and for physical states we will sum over the copies. A generic state of an 
electric charge-anti charge pair at vertices (v, v') will then be represented as a matrix 

|M R ;(v,vO) = ^£M« v (^= l ^ 1 |(^ p!e );(v J -))|(^, e );(v' J -))) 

y/\R\ M,V y V\ K \ P=0 ' 

with the normalization chosen so that L|fi=o l^^vl 2 = 1^1- Let's see how such a state arises in this spin lattice 
model. Recall that charge states are obtained by acting on a state j*?) at an intermediate stage of ground state 
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synthesis where the gauge symmetrization has been performed at all vertices except at v,- ;,v,- Applying the 
following joint charge projection operators onto j 1 ?) gives 



p« P (vij)p v %(vi J+ i)\v)/V¥Y\ 

ZhMeG[m%M h %3 T h^ij) T h<nj+i)\^}/V\\' 



R\- 
R 



L/iggLy 



^eG^(^ 1 )1 ft r^(^0lY,p^W]v,p^-i(vy^v ! ; ; - + i)|W)/y|H (22) 

r Lw-l6G^^(^ 1 )1ftAv^'-l(v«J,V/, / +l)|^)/^Fl 

$ Lw-ieGmh'- 1 )*)^,-! (v u , vy+i) W / vIR 

The operator 7/y^-i(vij,Vj j+i) = 7/i(v/ j)7V(v/ j+i) acts as LhR^r-iieu-jj+i) on the connecting edge. Before the 
projections, the system particle eij ; u+\ on the edge connecting the vertices is in state \e), while after the local 
gauge transformations on the boundaries it is in the state \hh'~ l ). Hence the charge state Eq. 21 can equivalently 
be obtained by beginning in the fully gauge invariant ground state \GS) and applying a projection operator that 
depends on the local state of the connecting edge e I j; , j+i: 

(23) 



l(^p>e);(vu ,-))\(^,e);(v iJ+1 ,-)> = VW\ L ls> Wl <SPterW \GS). 

g EG 

This argument in fact extends to electric charge anti-charge pairs separated by longer chains. Consider the chain 
[vi, V2, . . . Vyfc] and an initial state |\P) which is gauge symmetrized over all vertices but those in the chain . Applying 
the projections onto the charge pair at the boundaries and onto charge zero for the vertices {v2,V3, . . . \'k-i}, we 
have 

^(v< P (vi)m B^E C ! R (*i*j 1 )V^(vi I n)iA l ,„ Vl nS^ iJ (v;)|¥) 



The gauge transformations act on a state for the chain as 

^ lA -i(vi ! Vi)Eft 2) ... % _ i ny;2^(V;0Kl}[v I ,v 2 ]|4)[v 2 ,v3]---|4-l}[v t _ 1 ,v t ] 
= E*2,-A*-l \ h l^ h 2 )[vi,v 2 ]l^2^ 1 )[v 2 ,v 3 ]---|%-l4-l^ )[v jfc _ 1 ,v t ]- 

It is then the product of the group elements £\l2 . . . Ik- 1 along the edges of the chain that is an invariant under local 
gauge transformations on the vertices (excluding the boundaries). Notice that if one of the edges had opposite 
orientation to the chain, say e = [vy,vy_i], then the invariant would be l\ . . .£j-i£j l £j+i . . .£k-i- Hence we find 



that an arbitrary state of an electric charge pair along the chain [vi 

|(^p, e );(vi ,-)>l«p,e);(v* -)) = 



R 



■A-lGG 



,\'k] spanned by the edges {ef\ is 

k-i 



)Z\\hj)ej(hj\[R(UhrrU \GS), 

7=1 



(24) 



where o(e } ) = 1 if ej = [vy,v y -+i] and o(e } ) = -\ if e/ = [v/ + i,v y ]. 

We now describe how to construct the charge state of Eq. 21. This state does not violate any face constraints 
(since the local gauge transformations do not change magnetic flux), but it does violate the vertex constraints 
A(v,j),A(vy+i). To create this state, first prepare the vertex ancilla Vjj in state \e) Vjj . Apply the conditional 
unitary K(vij,ejj-jj + \) defined by 



K(v,e) 



L g e G \g)e(g\®R g (v) 

Z g eG\g)e(g\®R g -^>) 

and measure ancilla v,-j in the basis 



= [v, *J 
e = [*,v] 



\R 



ft,V/ 




Zmgywg) 



g eG 



The orientation of the edge dictates the left or right action taken by the operator K on the vertex ancilla. 

The above method only works probabilistically with probability 1/|G|. It is possible to prepare the vacuum 
electric charge states with unit probability as we know describe. For Eq. 19, the vacuum electric charge states at 
the boundaries of the edge ejj-jj + \ are 



l|*|;(v U ,v iy+1 )) = i£|(^ p , e );(v u -))|(^, e );(v u+1 -)). 
|K| ftp 



(25) 
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Using Eq. 23 this can be rewritten as the effect of a single qudit edge operator acting on the ground state, 

|l|R|;(vy,V ( y+i)) = L^LgeG^e^ij+AsimsTWlGS) 

= W R (e lJ;iJ+l )\GS) KM) 

where 

W R (e)=Z,\g) e {g\r R (8)- 

g£G 

These operators obey the sum rule, Y,r \R\Wr = \G\ \e) (e\. For S3 we have explicitly, 

w Rt = 1 6 

W R - = | e )( e | + |c + )(c + | + | c _)( c _|-| fo )(ro|-ki)(r 1 |-k2)(f2| (27) 
W Rl - 2\e){e\-\c+){c + \-\c-){c-\ 

Notice, that for the one dimensional irreps, the electric charge creation operator is a diagonal unitary acting on the 
edge. The operator W Rl is not unitary but we can construct it using adaptive measurements. First prepare the vertex 
ancilla vy in the state |+) Vj ; where |±) = (|0) ± |l))\/2 and apply the controlled unitary operation |0) v ,- ; -(0| <g> 
t/o (e^y;,-,.^ 1 ) + 1 1 ) v,,,- ( 1 1 «) C/i ) where C/ — diag( 1 , 1 , ^, ^* , ^, ^* ) , C/i = diag( 1 , — 1 , x / 2 , ^ 1 /2 ^* ^) in the 

basis {|e), \to), |fi),|?2)j \ c +)i l c -)}- Next measure the ancilla in the basis |±) v . .. The outcomes ±1 occur with 
probability p± = jT \ - For the outcome +1, the operator W Ro is successfully applied and the vacuum charge 
state is created. Otherwise we must correct. One strategy is to first apply the operator U%{e.i where U2 = 

diag(l, — i). This corresponds to a combined action on the system diag(0, 1,1, 1, — ^r,— ^)|GS). Now, 

re-prepare the ancilla in the state \e) Vij and entangle the edge with it by the operation X = £/=o \tj) ei - +1 (tj\ <8> 
Rtj ( v i,j) + 13 <8> 16- Next, apply a controlled gauge transformation at vertex vy that depends on the state of the edge 
via the operator Y = Lj=o |*/)ej h -+i (Ol ® T t -i (vy) + 13 <E> 16- This acts to map the three components of the code 
wavefunction with edge states in |f ; ) to the state \e). Since all these operations preserve the zero flux condition 
on neighboring faces, the states of the remaining lattice spins are disentangled from the \tj) components of the 
ancilla. Finally, we need to disentangle the ancilla from the lattice spins. This can be realized by first applying 
the single qudit operation on the ancilla U 3 = ( | e) ((t \ + (h \ + (?2 1 ) + |/i ) ( (to \ + \ + £,* {h \ ) \/3 + 1/2) ( {to \ + 
^*(fi I + ^(?2|)v / 3 + |fo) (e\) © I2 and then measuring the ancilla in the Fourier basis \k) Vi .. For outcome k apply 
the correction gate Z (eu-jj+i)- This then realizes W R2 (eij-jj + i)\GS). 

One can prepare states with electric charge particle anti-particle pairs that are separated further than one 
edge, e.g.|(P^ ,e); (v;, 7 -, — )) |(P*Q,e); (v, J+ 2, — )) as follows. Prepare vertex ancilla v;j in state \e) Vjj , apply the 
conditional unitary K(vij,eij-ij+i). Swap the qudits at locations vy and eu-jj+i and then swap qudits at locations 
e i,jj,j+i and Vjj+i. Apply the conditional unitary K(vjj + i,eij + \-ij + 2) and measure v,j in the basis {|^, v )}- To 
deterministically prepare the vacuum state \l\g\', (v,y, Vij+2)), apply the same sequence of swaps but instead of 
measuring the ancilla, apply the gate W R to the ancilla and then invert the sequence of conditional unitaries and 
swaps. This protocol easily extends to preparing arbitrarily long separated electric charge pairs provided care is 
taken during the controlled flip operations on the ancilla to keep track of the orientation of the edge for the system 
control. In the case of the non-unitary gate W Rr ,, another ancilla should be used to perform the operations on the 
ancilla which carries the information of the accumulated product of group elements along the path between the 
boundary vertices where the charges will be created. 

To move an electric charge from a vertex v to a neighboring vertex v' consider that the edges are oriented 
[*,v], [v,v']. Then we need to coherently map the state of the edge [*,v] to the product of the states of [*, v] and 
[v, v']. This is achieved by the following operator: 

A(v) E T g iv')T g (v)\ g } bi 4g\®\ g % v , ] ( g '\ 
g^eG 

The last projector A (v) reinforces the gauge symmetry at v so that total charge is conserved. More concretely, say 
we wish to move a electric charge one unit to the right: 

|(^p, e );(vy -)>-H(^ p ,«);(v u+ i -)) 
then the following protocol will work: 

• Prepare the vertex ancillae vy,vy+i in state \e) 
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• Coherently map the state at edge c,- to vertex ancilla v,- ; using the operator K (v; ,-,6; ,-) and 
similarly for the other edge using (v; ,+1,6; j ; / 

• Apply the controlled gauge transformations W(vjj) and W(v,- 

• Disentangle the ancilla v (J - by applying K(vu,eu-jj+i). 

• Prepare the vertex ancilla in state |0) v . .. 

• Apply the controlled gauge transformation W(vy). 

• Measure vertex ancilla in basis 

• For the outcome \J) apply the single qudit operation Z-ifa j-i-j j). 

• Disentangle the ancilla v,-j+i by measuring in the basis \Rfi,v)v lj+i ■ 

Moving charges in other directions is straightforward keeping in mind that if the edge orientations are reversed 
then the inverse gauge transformations should be applied. 

Finally, fusion of an electric charge pair \M R ; (v, v')) is realized by moving the charges until they overlap at 
one vertex (say v'). The outcome of measurement on the ancilla v' in the transport steps above determines the 
residual charge. For outcome \R^) = |0) the charges are perfectly annihilated into the vacuum, otherwise there is 
some residual charge. 

3.3.3. Dyons in D(5a) We briefly mention how to create dyonic excitations without discussing how to move 
them. Recall from Eq. 5, that there are three dyonic particles in D(S 3 ). The explicit representations of the dyonic 
particle antiparticle pairs located at {vij,fij) and (vjj+^/j-ij) are: 

^Ep^K^Up);^,//,;))!^ 1 ,^ 1 );^!,/;-!,;)), 

^Ip^K^.CpV.Cvy./y^KP^.Cp^iCvy+iJj-ij)), (28) 

^E}=ol(^.^);(vu./u)>l(^*.^);(vy+i,//-u)>. 

To create these states first create the vacuum magnetic charge states |0^j; (v;j,/y), (vy+ij/i-ij)) with [c] for the 
first two and [t] for the third state. Next apply one of the following projection operators 

W R{ = \(\e)(e\+!~\c + )(c + \+!~*\c-)(c_\), 

W Rj = W* { (29) 

w R 4 = £(k)<«|-|fo>fol)- 

These projections cannot be done unitarily with unit probability but similar to the implementations of Wr 2 one can 
use an adaptive protocol to realize them with unit probability. 

3.4. Anyonic interferometry 

Let's see how our spin lattice model reproduces these braiding relations given in Sec. 3.3. First consider braiding 
of fluxes. In Sec. 3.3.1, we described how to create a total magnetic charge zero state like Eq. 16. One process we 
could attempt is to simulate is to create two pairs of total charge zero states, braid a member of one pair around a 
member of the other, and annihilate. 
Say we begin in the state: 

J_ v \(r>Rt W„. . f. .\\\(vRt 



(30) 



l^n) = 75l P = ± l(^' ,c P );(v 2 ,4,/,, 3 ))|(/^ ,c_ p );(v 2 ,4,/ M )) 

^^=0,1,2 (V3,3,y 2j z))|(^,fi); (V3 j3! / 2l 4)> 

If we braid the flux [c] at /13 around the flux [t] at /2,2, we obtain the output state 

l^om) = ^E;=o,i,2lp= ± l(^ ,c_p);(v 24 ,/i,3)>|(P<,c_ p );(v 2 ,4,/ L 4)> 

using the relations: c p f ; c_p = tj+p and c p tjCptjC- p — Cptjtj+p — c p c £j . = c_ p . Now we can attempt to 
fuse the [c] flux anti-flux pair by propagating the flux at (vi,3,/i,3) one unit to the right. This is done by 
preparing the ancilla /i, 3 in state \e), mapping the flux to /i^ by A(v2.4,/i, 3 ) and applying the controlled unitary 
T(/i,3,ei,4 ; 2,4). If no braiding had been done by the [t] flux then the ancilla would be left in the state |(j) + ) f x 3 
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where 1^) = (\c+)±\c-))/y/± What we find is that the ancilla will not be disentangled from the system because 
the fusion is incomplete. Rather the ancilla is left in the mixed state p(/i,3) = j(\c + ) f l3 (c + \ + |c_)/ 13 (e_|). 
To measure this incomplete fusion, we measure p = (0~|p(/i.3)|(|) _ )- Before braiding, p = while afterward 
p = l/2. 

Next we examine the braiding of charges. Braiding of charges around each other acts trivially in the spin 
lattice quantum double model here because all the operations for creating and moving electric charges are diagonal 
in the logical basis and hence must commute. To obtain the other relations first consider the effect of applying a 
local gauge transformation to a electric charge anti-charge pair separated by a chain of vertices [vi , V2, • • • v k ] . For 
simplicity we consider that all the edges have the same orientation as the chain, i.e. e\ = \y\ , V2\ : e2 = [v2, V3] etc., 
though this is not necessary as in the case of an edge with opposite orientation to the chain one simply takes the 
inverse group element in the projection. The action is 



T h (vi) 



P^,e)-(v h -))\(^e)-,(v k -)) 

" l Lh lt ...h k eG \h) ei ® • ■■®\h k ) ek (h k \[R(h • . .**)*] J V> (vi)|GS) 



Th (vi 

Lh 1 ,...h k eG \hhi) ei (hhi\ <E) • • • <g> \h k ) ek (h k \[R{h\ 

Y*h'..h k eGWl)e,{ti l \® 



■®\h k ) ek {h k \[R{h-% 



■®\h k )e k (h k \ZAR{h- y ) 



h k )*]fi,v 
■■h)%.y 



\GS) 
\GS) 
...h k )Y 



(32) 



LARih-'rUKP^eUv, -))\(P%,e);(v k ,-)) 
Umk.m^livi ,-))\(I** ,e);( Vk ,-)) 



\GS) 



where we used the gauge invariance of the ground state. Similarly we can compute the action of a gauge 
transformation on the anti-charge, and we find for a general electric charge vacuum state: 

7),(v)|M*;(v/)) = \R(h)M R ;(v,V)), T h (v')\M R ;(v,v')) = \M R R{h~ l ); (v/)>. (33) 

But this local gauge transformation Th(v) can also be viewed as creating a flux anti-flux pair h,hr l and braiding 
it in a counterclockwise sence around the electric charge at v followed by annihilation of the flux pair. The 
annihilation probability for the flux was zero, meaning the only non trivial action was on the electric charge 
pair as indicated, hence our model reproduces the correction braiding relations for a flux around a charge. In 
particular, given that we can create the vacuum state for the representation R: \1\r\', (v,-j, v M )), we can create the 
state \R(h); (v,j,Vp i9 )), by applying 7),(v;j), which is a product of four single qudit unitaries. How can be measure 
the effect of this braiding? One way to is invert the vertex ancilla steps that prepared the electric charge pair. For 
the state |1iri; (v,v')), the inverted sequence will leave the ancilla disentangled from the system, whereas if the 
initial state is \R{h)\ (v,v')) then the ancilla will be left in a mixed state. 

To actually extract the value of amplitudes for fusion of charges into the vacuum we can use the vertex 
ancilla prepared in |\|/+) where = (|e) v ± |/z) v )/v / 2, and apply the controlled operation W(v) followed by 
measurement of the ancilla in the basis with outcome m = ±1. The outcome distribution satisfies 

9l[tr[R(A)] 



P(m=l)-P(m=-l) = X[{l lR r,(v,v')\R(hy,(v,v')}] 



\R\ 



which is the real part of the fusion amplitude for R(h) — > vacuum. Similarly, measuring the ancilla in the basis 
|\|/*) = (|e) y ± i\h) v )/V2, yields the imaginary part of the fusion amplitude for R(h) — > vacuum. 



3.5. Elementary operations for quantum computation with anyons 

We now show explicitly how to create anyonic states and perform braiding and fusion operations which are 
universal for computation. This section follows closely the work of Mochon [18] who proved two important 
facts: first that by working with magnetic charge anyons alone from nonsolvable, non nilpotent groups, universal 
quantum computation is possible, and second that for some groups that are solvable but not nilpotent, in particular 
53, universal quantum computation is also possible if one includes some operations using electric charges. 

The first step is to identify the logical basis. We will encode a qutrit in the three charge [t] magnetic fluxes: 
to,ti,t2- We will work with pairs of flux with total trivial flux and will adopt the simplified notation for pure 
magnetic charges located at positions r and s to \tj,tj ,(r,s)) = \(P Rl ,tj);(v,fj))\(P R i ,tj );(y,fj)) and for an 
electric charge pair at positions r,s, \M R ;(r,s)) = \M R ; (v r , v s )). Since we will be working with magnetic fluxes 
that may not have zero charge, the base point vertex v is important but for convenience we suppress it keeping in 
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mind that all fluxes should be valued with respect to some fiducial base point. The Z basis is {\tj,tj ')} and the X 
basis is {\tj,tj~~ )} where 

\f J J j - l -(r,s)) = -L^\t k ,t^-(r,s)) 
v 3 

The state |fo,fo _ ) is just me vacuum magnetic charge state |0[ ( j). The generators of the Pauli group operations on 
this basis are 

X(r,s) = t lO+i,*^; (r,*)> (tj,*/ 1 ; (r,a)|, Z(r,s) = £ %%tj\ {r, S ))( tj ,tj\ M| 

7=0 7=0 

Encoded information is always stored in flux anti-flux pairs, and the braiding operations we employ will always 
braid the pair together thus perserving the total zero flux condition (see. Eq. 3.3) though there may be a residual 
electric charge. First we describe how to initialize single qudits in the X or Z basis. The initialization of the 
state \tj,tj~ ) is equivalent to preparation of the vacuum state described in section 3.3.1 except where we alter 
the correction step via the face ancilla / accordingly, i.e. we apply [1|g|-|[£']| ©Zp^K/) given a measurement 
outcome \k^)f. To initialize the state \tj,tj , {f\ifi)) we require a means to perform a projection. The main 

ingredient necessary is the projection of the qutrits onto the subspace ac > = span c {\tk,t k ),k ^ j}. This can 
be accomplished with the assistance of an ancilla pair in the vacuum two dimensional electric charge state 
|1r 2 > ( r i > r 2))- Say the initial state of the qutrit is 

1^(1,2)) = £ c ^,/7 1 ;(i,2)) 

7=0 

and we want to project out the to component, i.e. we want to project onto K^'o . To do so we braid flux 1 around 
the ancillary electric charge at r 2 creating the correlated state 

2 

^ r2 \^)\l R2 ;(r u r 2 )) =^cj\tj,tr l ;(l,2))\R2(tj);(rur 2 )) 

7=0 

Next we apply the ancilla assisted local gauge transformation T f -t (r 2 ) so that the joint state is 

c |fo,f -^(l,2))|l« 2 ;(r I ,r 2 ))+ Cl |f 1 ,ff^(l,2))|/? 2 ( C+ );(r 1 ,r 2 ))+ C2 |f2,^^(l,2))|/? 2 ( C _);(r 1 ,r 2 ^ 

Then we annihilate the electric charge pair and observe whether it fused into the vacuum (using the vertex ancilla 
protocol described in Sec. 3.4). The states |/? 2 (c±)) always fuse into one dimensional representations (either 
the vacuum state or the sign representation each with probability 1 /2). If the charges fuse into the vacuum 
then we throw the qutrit out, whereas if the outcome is imperfect annihilation then the projection onto the 

subspace is successful. Similarly, to project onto oc'j we use the same protocol but with the local gauge 
transformation T r -i(r 2 ). By composition of projections we can prepare any basis state beginning in |0[,j), e.g. 

|»i,»r 1 ;(l,2))=i:^3C^|0 w ;(l > 2)>/vlHl- 

Let's see how to generate the controlled sum operation £((1,2); (3,4)), where 

2 

E((r,j);(p >? )) = £ \t h t-T x -(r,s)){t h t^-{r, S )\®X>( PA ) 

7=0 

To realize this we need two primitive braiding operations. First consider the anicilla assisted operation with the 
ancilla prepared in | ?o , f/T 1 ; (5 , 6) ) , 

o((3,4);(5,6))[^ 2 3^ 2 4 ]K 3 2 5 ^ 2 6 ]|f i ., f ,- 1 ;(3,4))|ro,^ 1 ;(5,6)) 
= o((3 , 4); (5, 6)) [% 2 3 ^ 5 2 4 ] \tk, h 1 ; ( 3 - 4))^ 2i |f , 1 ; (5 , 6)) 

= o((3,4);(5,6))|%0% 1 ' f 2A 1 % 1 ;(3,4))|f 2 ,,f 2 - 1 ;(5,6)) (34) 
= |f 2 ,,f 2 V;(3,4))|f , fo - I ;(5,6)) 

The final swap o on the two pairs is just a classical operation since each pair has total charge zero. Second, 
consider the action of braiding one control flux around the target pair 

^ 3 ^ 2 4l^?7 1 ;(l ! 2))|^,%- 1 ;(3 ) 4)) - l^r 1 ; (1,2)) \t2j+2k,t£ +2k ; (3,4)) (35) 
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The composition of these two operations realizes the inverse controlled-sum gate: 

IO,0 1 ;(l,2))|r,, fA 7 1 ;(3,4)) -> \ t] ,tf\ (1,2)) \t 2j+ 4k,t 2 ' ; (3,4)) 

= |^,/7 1 ;(l,2))|/ 2y - +fc) f-^;(3,4)) (36) 
= [E((l,2);(3,4))]- 1 |o,0 1 ;(l,2))| f ,, f ,- 1 ;(3,4)) 

Now [E((1,2);(3,4))- 1 ] 2 = E((1,2); (3,4)), hence by two applications we obtain the controlled sum gate. Using 
a controlled-sum gate with the control prepared in state \tj,tj ; (1,2)) allows the implementation of an X-> gate 

on the target. Imagine that we had an ancillary source prepared in \fi,fj~ \ (3,4)). Using a controlled-sum gate 
with the target being this ancilla we can perform the one qudit operation Z 7 on the control which is verified by 
considering the action on a complete basis for the control: 

E((l,2);(3,4))| f l,r 1 ;(l,2))|f^fr^(3,4)) = [Z^ f l,r^(l,2))]|0,0- I ;(3,4)). 

How do we prepare such an ancilla? We can do this by preparing two ancillary pairs in the state 
\ x ¥) = |0[,j; (1,2)) \to^Q l ; (3,4)) and applying the controlled sum gate £((3,4); (1,2)). This prepares a 
maximally correlated state and if we trace over the pair (1,2), we obtain the maximally mixed state P3.4 = 
i^ = ol?y,?7 1 ;(3,4))^ ) ?7 1 ; (3,4)| = Il2 =0 | r ~, f -i;( 3i 4))( , f 7i;(3,4)|. Now using the methods of [18] one 

can measure pi 2 intheXbasis which will project the state onto |f~ ,f) _1 ; (3,4)) with probability 1/3. This process 
can be repeated until the desired outcome is obtained. Such a measurement requires only the preparation of many 
ancilla pairs prepared in |0[,j) Similary, one can measure in the Z basis with the assistance of many ancilla pairs 
prepared in 1 1 , 1 ) . 

Up to now we have shown how to implement the Clifford group operations: controlled-sum, and X a Z b as 
well as preparation in an X or Z eigenstate and measurement in the X and Z basis. In order to realize universal 
quantum comptuation we require the ability to perform a non Clifford operation such as the three qutrit Toffoli 
gate defined by the action: 

r((l,2);(3,4);(5,6))|^*7 1 ;(l,2))|/ k ,/r- 1 ;(3,4))|//,^ 1 ;(5,6)) 
= |/,-,/7 1 ;(l,2))|/ fc) ^ 1 ;(3,4))|^ + ,,^;(5 ) 6)) 

Because the group S3 is solvable we cannot implement the Toffoli gate by braiding magnetic charges alone [27]. 
However, as shown by Mochon [18], we can use allowed braiding operations with electric charges to aid in 
preparing magic states which facilitate the Toffoli gate. We require two types: 

I <>Afi ; (n , n) ; {n , n); (r 5 , r 6 )) = ± E^=o \tj,tj l \{n,r 2 ))\t k ,t- l ;{n,r A ))\ t jk , tj k l ; (r 5 , r 6 ) ) 
ftm; (n .rz); (n, *)> = 3 E?jmj £ 5 °' Ao l'y>*7 u > (i> rz)) Kh u > , r 4 )) 

To realize Toffoli gate on the intial state of three charge pairs 

2 

|¥;(1,2;3,4;5,6)>= £ %, b , c \t a ,t- 1 ; (l,2))\t b ,h\ M)^ 1 ; (5,6)) 

a,b,c=0 

we append the magic state \§mi', (n ,r 2 ); (r3,r 4 ); (rs,r^)) and apply 

[£((5,6);(r 5 ,r 6 ))][£((r3,r 4 );(3,4))]- 1 [£((n,r 2 );(l,2))]- 1 

followed by measurement of 1 and 3 in the Z basis with outcomes 1111,1113 and 5 in the X basis with outcome 1115 
and subsequent correction gates X(r\ 1 r2) mi X{rj,,rn) mi X(rs 1 r^)^ mi ' ni . The outcome on the ancilla states is 

2 

£ PaAc^K'C 1 ; (n,r2)) \tb,h l ; (r 3 ,r 4 )) \t a b-m 1 b-m 3 a+c,t~ b 1 _ mib _ m3a+c ; (r 5 ,r 6 )) 

Finally we apply [Z(r 5 ,r 6 )]-'"5 [£((r 3 ,r 4 ); (r 5 ,r 6 ))p [£((n,r 2 ); (r 5 ,r 6 ))] m 3 to realize 

|X> = t P^r" ,5a6 k a ,^^(n,r 2 ))|f & ,f i -^(r3,r 4 ))|U +t -,C c ;('-5 I '-6)) 

This set of operations has the action of the Toffoli gate acting on up to a phase ^ m ^ ab on the first two qudits. 

The correction gate C = Y^ab^^a-,^ ki>>%~ ) ('a,*" 1 1 I) i s m tne Clifford group and one might 

hope that it could be generated from our set of available operations, e.g. using the techniques in 
[28]. However, to do so requires the single qudit discrete Fourier transform gate which is not accessible 
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using braiding operations. Hence we use the second magic state \§M2'> 0*7: r 8; r 9) r io)) an d a PPly the gates 
[E((r3,r4);(r9,rio))]" 71 P((ri,r2);(r7,r8))]'" 3 followed by measurement of rj and rg intheZbasis with outcomes 
niT,m<). The outcome is 

a,b,c=Q 

We can append another magic state \§M2\ ('"H,'*i2; r i3)'*i4)) and measure again with outcomes mn,/ni3. Repeating 
k times the total phase accumulated is 8(k,a,b) — Y!x=\ ^ m 3+4/b a ^>"5+4i> ^ n or der to realize the Toffoli gate we 
demand that for some k, 8(k,a,b) = m^ab mod dVa,b, where d is the qudit dimension (here d = 3). Each step of 
appending a magic state, applying controlled operations, and measuring applies a phase on the state. This is akin to 
picking k balls with d 2 "colors" which are labeled by the pair (a,b) from an independent and identically distributed 
probability distribution and placing then into d 2 bins each labeled by (a,b). We find a satisfying distribution if the 
number of balls in each bin n a j, = m$ mod d\/a,b. The probability that after k trials the condition is not satisfied 

is (for d prime) approximately (1 — l/d cl2 ) k . The preparation of the magic states follows identically the procedure 
in [ 1 8] with frequent use of the projection onto 1Q *o . 

4. Physical Implementation in atomic spin lattices 

An experimental candidate for realization of the ideas herein is cold trapped atoms in a two dimensional optical 
lattice (with motion quenched in the third dimension). Single lattice site occupancy of atoms prepared in motional 
ground states can be prepared in a region of an optical lattice using a variety of loading techniques [29]. In 
our construction we use an antiferromagnetic array of two type of particles: system particles or A species, and 
ancillary particles or B species which will assist in state preparation and generation and dynamical propagation of 
anyonic excitations. This bipartite division could be simply a spatial labeling of the same type of physical particle 
or it could correspond to physically different atomic species or subspaces of the same species. The advantage 
of using different species for system and ancilla particles is that it could assist in addressability of the controlled 
unitary operations and measurements in our simulation. In Fig. lb we suggest a suitable 2D lattice architecture for 
trapping and manipulating the atoms. By tuning the intensities and phase of the trapping fields, particles can be 
brought together pairwise with left, right, top, and bottom neighbors to facilitate controlled gate operations. Single 
particle coherent control and pairwise entanglement generated by controlled collisional exchange interactions have 
been experimentally demonstrated with this architecture [30]. 

The six group elements of S3 can be encoded into the ground electronic hyperfine states of a trapped alkali 
atom with enough available levels. Candidates include 87 Rb or 23 Na each having 8 available ground hyperfine 
states or 133 Cs having 16 levels. Arbitrary single qudit operations can be done using Raman laser pulses or 
microwave pulses that connect pairs of hyperfine states and measurement can be done using a cycling transition 
(see Fig. 2). The atoms are prepared in an optical lattice with spacing X/2 where A, is an optical wavelength. At 
this spacing and for tight lattices, the atoms are well localized and do not interact. In order to perform controlled 
unitary operations, one strategy is to bring atoms closer together pairwise and turn on a magnetic field which 
would generate a trap induced shape resonance on a pair of internal states. Ideally one would like this collision 
to be stable in the sense that it would be diagonal: H mt = U\gj)A (gj | <£> \gf) b (gf | ■ This could be provided for by 
choosing the collision on maximal angular momentum projection states of the ground hyperfine level, chosing e.g. 
\gj) = \gf) = \F^,niF = — Fj.) since a symmetry of the dominate collisional interaction is mp A + mp B . Quantum 
gates using such diagonal collisional interactions generated via trap induced state resonances have been proposed 
in [32, 33]. There it was shown that a robust controlled phase gate between a pair (A,B) of 133 Cs atoms is 
generated by collisions between the maximally projected angular momentum pair \Ft,mp = Ft^lFi ,mp = F^)b at 
a particular lattice well spacing. In order to have high fidelity gates, it is important to work in a parameter regime 
where there is negligible interactions between other states for two reasons. First, non maximally projected angular 
momentum states have elastic collisions which are non diagonal interactions, making a target two qudit unitary 
difficult to achieve, and second inelastic collisions between states in the upper hyperfine manifold Ff can destroy 
internal state coherences and convert the internal energy into kinetic energy of the dimer. 

Another possbility is to use the collision gate suggested in [34]. This would work by first using Raman 
pulses to map the states \gi)AB to the first excited vibrational state of lattice wells A and B. Then adjusting the 
intensity of the trapping fields such that the well minima are tilted with respect to each other and bringing the wells 
together (but not overlapping) pairwise and finally removing the bias in the wells. At an appropriate well minima 
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Figure 2. Proposed implementation used atoms trapped in a optical lattice, (a) In the 2D retroreflected optical lattice 
of Ref. [38], when the trapping laser field polarization is in plane, then the lattice has periodicity X/2 whereas if the 
polarization is out of plane, the periodicity is X with a relative shift of X/2 resulting in tunable an array of double well 
potentials. The upper panel shows the optical lattice for zero field component with out of plane polarization. In the lower 
panel, nearest neighbors along one dimension are brought together pairwise by increasing the intensity of the field with 
out of plane polarization. Left neighbors or right neighbors are joined depending on the sign of phase of the out of plane 
polarization light. Top and bottom neighboring pairs can be joined by increasing the out of plane polarization intensity 
and changing the phase of the field component with polarization in plane. We note that this architecture allows for global 
addressing of the A or B species by virtue of distinguishability in the motional states of the two type of lattice wells, (b) 
Encoding a d = 6 qudit into the ground hyperfine levels of a 87 Rb atom. The logical states are labeled by group elements 
of S3. In the figure, the states |c+) and fa) are shown coupled by a Raman laser pulse with Rabi frequencies 12 1 2 by the 
Hamiltonian Hal = ^eff(^|c+){f2 + e~'^\h)(c+\), where fl c ff = J2iI22|/A and (j) = argfiliilj). Arbitrary unitaries 
in the two dimensional subspace can be generated from this family of Hamitonians, and given a connected coupling 
graph between basis states, which is provided for here using polarization and frequency selectivity of Raman pulses, 
any U 6 SU(d) can be generated [23]. Projective measurements can be made by mapping population in each sublevel 
to \Ff = 2,mp = 2) and measuring resonance scattering on the cycling transition \F> = 2,mp = 2) — > \F' = 3,mp = 3) 
transition. The state |c_) can be used for entangling gates via ground state collisions, described by the interaction 
H mt = U\c-)a(c- \ ® |c_)b{c_| between neighboring atoms A and B. This interaction is diagonal due to total angular 
momentum and energy conservation. The state \r) is used as an ancillary state to shelve amplitude in the state |c_) of 
atoms that we do not want to participate in a controlled interaction. 



separation, there will be a large tunneling and large collisional shift acquired on atoms in the first vibrational states 
which are nearly degenerate and near zero tunneling shift for ground vibrational state atoms which are shifted in 
energy by an amount large compared to the collisional energy. Then the wells can moved apart again and the 
population in \gj)A,B mapped back to the vibrational ground state. Whatever entangling gate is used, given the 
ability to generate a phase on one basis state pair, arbitrary two-qudit gates can be generated [23]. 

In order to perform the quantum simulation here, some amount of addressability of the atoms will be 
necessary. This is slightly non trivial as the atoms are separated by X/2 and therefore cannot be directly addressed 
by optical fields. One solution is to use gradient field spectroscopy wherein a strong local electric field gradient 
(created e.g. using optical tweezers) is applied in the vicinity of the target atom. This field will also have some 
crosstalk with neighbors of the target atom but will shift the energy levels strong enough at the target location such 
that chirped microwave or Raman operations coupling hyperfine states act trivally at any location other than the 
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target [35]. To implement two qudit gates on a target pair of atoms A and B we suggest the following strategy: 
first for all atoms except A or except B map population in state \gj) to the spectator internal state |r), second 
tune the lattice trapping fields to bring A and B together and apply the unitary U = e '^ g j*> A ( g j^ 8 j*i B ( s j\ in parallel 
on all nearest neighbor pairs, third invert the mapping to return population to \gj). The invertible mapping on 
non target atoms can be done using STIRAP pulses with shaped intensity profiles as proposed in [36, 37]. For 
example, using a pair of Raman fields that couple states |c_) — * \r) (in the encoding of Fig. 2) one could make the 
Raman pair out of standing waves of light with a node at the target atom A or B. For any atom not at a node, there 
will be an invertible adiabatic mapping |c_) — > \r) which works with fidelity approaching unity. This technique is 
particularly well suited to addressing a line of atoms using a ID standing wave for the STIRAP pulse. One could 
also use gradient field spectroscopy to perform the mapping although this may not be as high fidelity due to the 
cross talk with neighboring atoms. 

5. Conclusions 

We have shown how universal quantum computation by anyonic braiding can be demonstrated in an optical 
lattice using the method introduced in [17], based on the creation and manipulation of anyons via controlled 
interactions of a code lattice with an ancillary species. In addition, we show how to construct the relevant 
topological state using the same method in absence of the suitable Hamiltonian. The system considered is the 
D(Si) quantum double lattice model introduced by Kitaev [19]. In all cases, protocols are thoroughly described 
and the experimental techniques required are discussed in detail. 

The experimental realization of this scheme relies in methods already in the literature. For the procedures 
described, which exhibit optimal scaling, the auxiliary species is distributed in a lattice with the same spacing as 
the code lattice, so addressability is an issue, but we propose ways to overcome this problem. We believe that these 
techniques are of sufficient interest to warrant experimental work. In particular, non-Abelian anyonic statistics can 
be demonstrated with our interferometry protocols. 
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Appendix A. The quantum double of a finite group 

Given a finite group G = {gi}, the group algebra C[G] is a complex vector space spanned by basis vectors gj 
(which can be chosen orthonormal to define an inner product.) Indeed, the multiplication gi ® gj i— > gjgj inherited 
from G makes C[G] an algebra. The unit element e can be viewed as embedding the complex numbers in C[G], via 
X i ► Xe. We can also define the algebra J (G) of complex functions on G, which has a basis given by Kronecker 
deltas, P/, : 8/ Ig . The multiplication in T (G) can be recovered from PhPy = 8/, yPj,. The unit is the map 
1 : g i-> Vig. 

Both C[G] and J (G) can be endowed with Hopf algebra structure. A Hopf algebra A is an algebra where: 
(i) Apart from the associative multiplication m : A (g>A — ► A and unit r| : C — ► A just described, there exists a 
comultiplication A : A — > A ® A and a counit e : A — > C (a coalgebra structure) satisfying (A <E) id) A = (id ® A) A, 
and (e®id)A = (id® e)A = id, where id is the identity map on A. For C[G], these maps are given by 

Mg)=8®8, E(g) = l, (A.l) 

while for f (G) 

A{P h )= £ Ph,®Ph 2 , £(P h )=8 h ,e- (A.2) 

h\h2=h 

For a general element, 

A(Atf)= L P^S^PhJ, E(P h g) = hA- (A3) 

h\hi~h 
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(ii) There exists an antilinear mapping S : A — > A, the antipode, satisfying m(S <E> id)A = m(id g) 5)A = r|e. The 
antipodes of C[G] and f (G) are determined, respectively, by S(g) = g~ l and S(Pf,) = P/ 



hr 



Moreover, the structures of algebra and coalgebra of C[G] and J (G) are dual to each other. This allows to 
use Drinfel'd's quantum double construction to define a quasitriangular Hopf algebra structure on J (G) x C[G] 
with multiplication Phg®Ph'g' ^ §h,gh'g- lP hgg' unit le, comultiplication P h g h-> Y*h x h z =hPhig ® Ph 2 g> counit 

i — > 8/ I ( ,, and antipode i— » /".-l^-i-g -1 . This is precisely the quantum double D(G) of the finite group G. It 
is generated as an algebra by the elements {lg\ g £ G} which form the electric gauge group G, and {P g e\ g £ G}, 
which are projections onto the set of states with flux g in the theory. 

The representation theory of D(G) is well known. Each irreducible representation is determined by a 
conjugacy class [£] = {ulu~ x \u £ G} in G, and an irreducible representation R of the centralizer of an element 
i £ [£], Nm = {u £ G\ ul = £u}. Let us label the different elements of the conjugacy class [a] as 



[a} = {hW,hW,...,hW} 

[a] 

The centralizer Ni a j C G, will be defined as the centralizer for the first element h\ . We can relate the different 

feci 

elements of [a] in terms of the first element h\ , via a set of representatives of the equivalence classes of G/N]^ 
y[a] _ r M [a] [a] , , [a] _ [a] , [a], [a] 

~ a'™' a'™' ' ' " a'™' _ a!™ 1 1 ^ a!" 1 ' 
and adopt the convention that x ^ = e. The carrier space for the irrep ^ is 



^ [a]) = {|Ar J ^)|Ar J e[a],0<;<|/?|-l}. 
The action of the irrep of an element of the Hopf algebra on this space is 

U m a] )( P Hi)\ht\vj) = \ gh n g -X\sht ] g-\R{g) m j v£>. (A.4) 

The element g = (x' a L )~ l gx l '*L is constructed from the gauge transformation g and the flux hf^ to satisfy 

g h\ > g - 1 h\ 

[g,h\ a '] = (verified using the fact that (x j , g -1 * j , = ftj ) implying that g £ N a y In this way 

feci 

we see that the the action of the irrep ^{Ph ® g) on the carrier space is to perform a gauge transformation on 

the centralizer charge followed by a projection onto the flux conjugated by g. The dimension of the carrier space, 

[ex] 

also known as the quantum dimension of the particle j is 

The quantum dimension satisfies the sum rule 

L(4Mj 2 =Ll[«]| 2 Li^(%])i 2 = Ei[«]| 2 i%]i = Li[ a ]|i G i = i G i 2 - 

[a] R a a 

There is some arbitrariness in the choice for representative for the equivalences classes and the ordering of the 
elements therein, however different choices lead to unitarily equivalent representations of the quantum double. 

The operation of braiding two particles is generated by the monodromy 3^, which effects a counterclockwise 
interchange of two particles: 

* = aoU m a] ) ® u %w)(L G p s®Phg)- 

Here a is the particle interchange operator and the operator in parentheses is the universal /^-matrix, an element of 
D(G) x D(G) describing a gauge transformation on the second particle by the flux of the first. The explicit action 

on the two particle state space ^ ®^r'(jv , ) i s: 
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where h\ 1 = (x , , „ , , )~ h\ x l , 1, is defined as above. Frequently, we are interested in the action of a full 

counterclockwise braid, of one particle around another. For the case of a pure flux braiding another pure flux: 

^ 2 1 /z , 0) | , 0> = | /z, [ 3 ?' ] )/z ^ /,[?' ] ) - 1 , 0) |/z I 011 /zi?' ] (/z ^ ) - 1 , 0> 
and for a pure flux braiding a pure charge: 

^ 2 i/ 1 | a ],o)i e ,v^)=^i/ 1 | ci i,o>| e ^(/ J | oi i) <> vf). 

I 

Appendix B. Some properties of the group S3 

S3 is the group of permutations of three objects, that we label {0, 1,2}, and is the smallest non- Abelian group. 
Elements of S3 are organised in three conjugacy classes, namely: 

• Identity e. 

• Transpositions (reflections) to = (01), t\ = (12), to = (20). 

• 3-cycles (rotations) c + = (012), c_ = (021). 
The multiplication rules for S3 are as follows: 

Uti = e, tjt k = c £j k for j^k, 

U c p ~ ti+p> c ffi = h—pi 

CpCp = c_ p , c c c T ~e for a 7^ T. 

together with the trivial operations involving e. Here e,- 4 = ± is such that k = j + e ;J t (modulo 3, as indicated.) 
In particular tf = tj, c^ 1 = c_, and conjugation relations are 

tititi = ti 1 

tfytj = tj where all of i, j, k are different, 

c pU c -p = *('+P) 
HCpti = C_p 

CpCgC-p = c a . 

The group has three irreducible representations (irreps): the two one dimensional irreps are the trivial one 
R~l(g) = 1, the signature representation 

R;(e) = +l=R-(c p ), Rj(t t ) = -1, (B.l) 

and the two-dimensional irrep 

R 2 (e) = l 2 , R2(tk) =a r exp ^/y ka z ^j ,R 2 (c p )=e\p ^'Py(r'j , 



Explicitly, 

/ 1 n \ / V- n \ 

R 2 (c-) = 
, Ri(t 2 ) = 



£,* 
% 
% 
^ 

where \ = e'' 271 / 3 , t,* = e'' 4lc / 3 . 

The characters %R(g) = tr [/?(#)] all equal to ±1 for the one dimensional reps and 

%R 2 (e)=2, XR 2 (*/)=0, Xr 2 (c p ) = -1. (B.2) 
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The permutation representation for 63 is a set of 6 x 6 matrices that faithfully represents group left action 
on the basis ||e), |?o), |fi), \t2), \c+), |c_)}, i.e. L/,\g) = \hg). Similarly, we have unitaries for right multiplication 



SaIs) = \gh 


>. The unitary matrices satisfy [L/ n /?/,/] = 


= and are 


given by 
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Actually, g 1— * L g and g 1— > R 1 determine the left and right regular representations of G, respectively. This is 
because right multiplication inverts the order of group multiplication R g R g i = R g ' g whereas R g = R -1 defines a 
representation. 

Finally, we point out the group S3 has a semi-direct product structure which may be exploited to simplify 
physical realizations. Recall the definition of the semi-direct product. Suppose that we are given a group G 
with a normal subgroup N, a subgroup H, and the property that any g E G can be written g = nh for n€ N and 
h £ H. Let (j) be the homomorphism (j) : H — > Aut(iV) where 0/j(n) = hnh . Then G is isomorphic to the semi- 
direct products xui? and the isomorphism identifies the product nh G G with the pair (n,h) (zN We have 
S3 = Z3 x^li2 = (a,b\a 3 — e.b 2 = e,b q ab~ q = a 2q ). Here the homomorphism is specified by 06 (a) = bab~ l = a 2 . 
Using the notation above we can choose Z3 = {e,c + ,c_} and 1^ = {<?,fo}, and any element g G S3 can be written 
g = c' + tQ for r G {0, 1,2}, s G {0, 1}. Introducing the basis for group elements {|r)|i) = |cj_?q}}, i.e. a product 
basis for a qutrit and qubit, we have compact representation of the left and right action operators: 

L e = 1 3 ®1 2 , L t0 =F{\,2)®o x , Lt 1= F (0,2)®^, U 2 =F(0 ! l)®cr x , 

L c+ = X~ 1 ®1 2 , L c _=X<g>l 2 , R e = h®l 2> /?, = l3®o x , m « 

R tl = X- l ®0- +X®o+, R h =X- l ®G + +X®a-, ( ' 

R c+ = X® I0X0I+X- 1 ® |1)(1|, R c _ =X- l (g>\0)(0\+X(g>\l)(l\ 

where F(i,j) = + ffi 1 flips two basis states of the qutrit. The electric charge creation operators of 

Eq. 3.3.2 assume a particularly simple form: W R - = I3 <g> O z and Wr 2 — diag(2, — 1, — 1) <£> |0)(0|, as do the dyonic 

projectors: W R i = diag(i, |, ® |0> <0|, W R z = W* u W R 4 = \0)(0\®o z . 
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